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Abstract 

A complete numerical calculation of the temperature anisotropies and polarization of the cosmic 
microwave background (CMB) in the presence of a stochastic helical magnetic field is presented 
which includes the contributions due to scalar, vector and tensor modes. The correlation functions 
of the magnetic field contributions are calculated numerically including a Gaussian window func- 
tion to effectively cut off the magnetic field spectrum due to damping. Apart from parity-even 
correlations the helical nature of the magnetic field induces parity-odd correlations between the E- 
and B-mode of polarization (EE) as well as between temperature (T) and the polarization B-mode 
(TB). 



1 Introduction 

The existence of magnetic fields in the universe has long been established by observations on 
small upto very large scales, that is on the scales of galaxies, galaxy cluster and superclusters [1]. 
Recent observations even indicate evidence for magnetic fields at truly cosmological scales [2] . Over 
recent years there has been an increasing interest in cosmological magnetic fields, from generation 
mechanisms to observational imprints on the cosmic microwave background (CMB) (for recent 
reviews, e.g., [3]). 

For a primordial magnetic field to serve as an initial seed field to explain the galactic magnetic 
field, it is assumed to be generated either in the very early universe, such as due to the amplification 
of perturbations in the electromagnetic field during inflation (e.g. [4]), or after inflation has ended, 
e.g., during a phase transition such as the electroweak or QCD phase transitions in which case it 
has been shown that the resulting magnetic fields are helical [5]. The generation of helical magnetic 
fields during infiation has been studied in [6]. However, in any case it is natural to assume that the 
magnetic field existed long before initial conditions are set for the evolution of the perturbations 
determining the temperature anisotropies and polarization of the CMB. Thus observations of the 
CMB could in principle be used to put limits on the parameters of a primordial magnetic field. The 
aim here is to calculate the angular power spectra of the temperature anisotropies and polarization 
for a Gaussian stochastic helical magnetic field. In this case the magnetic field two-point correlation 
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contains a symmetric and an asymmetric part. Each of which are characterized by one amphtude 
and one spectral index. The hehcal part results in non-vanishing correlations of the temperature 
anisotropics and the polarization B-mode on the one hand and of the polarization E-mode and 
B-mode of the CMB on the other hand. These have been first considered in [7], where the vector 
modes were studied. The spectra of the CMB anisotropics due to the tensor mode have been 
treated in [8] and vector and tensor modes in [9]. Moreover, the helical part also contributes 
to all other correlation functions, that is the temperature and polarization E-mode and B-mode 
autocorrelation functions as well as the temperature polarization E-mode cross correlation of the 
CMB. Here a consistent numerical treatment is presented taking into account the magnetic field in 
the initial conditions and the evolution equations as well as calculating numerically the correlation 
functions encoding the contribution of the magnetic field. Using a Gaussian window function the 
magnetic field spectrum is effectively cut-off at a damping scale. Moreover, the angular power 
spectra of the temperature anisotropies and polarization of the CMB are calculated for scalar, 
vector and tensor modes. 

The calculation of the CMB temperature anisotropies and polarization requires to solve the 
evolution equations of the perturbation equation of the geometry and matter components together 
with the Boltzmann equations for photons and neutrinos. Since the first numerical program COS- 
MICS [10] to do this and with a significant improvement in speed using the line-of-sight integration 
in CMBFAST [11] there has been a steady evolution of numerical codes such as CAMB [12], CM- 
BEASY [13] and CLASS [14] being the latest addition. The effect of a primordial nonhelical mag- 
netic field on the CMB anisotropies has been calculated using different approaches: synchronous 
gauge and thus a modified version of CMBFAST in [15, 16] or the covariant formalism and a 
modified version of CAMB in [17, 18]. Finally using the gauge- invariant formalism for the scalar 
perturbations the CMB anisotropies have been calculated in the presence of a stochastic nonhelical 
magnetic field using CMBEASY in [19]. 

Here the modified version of CMBEASY [19] has been expanded to include firstly the numerical 
calculation of the correlation functions of a helical stochastic Gaussian magnetic field and secondly 
a new part to solve the Boltzmann equation and calculate the CMB temperature anisotropies and 
polarization for vector and tensor modes. These are calculated using the total angular momentum 
approach of Hu and White [20] . In section 2 the decomposition of the magnetic field contribution for 
scalar, vector and tensor modes is described. Section 3 is devoted to the calculation of the relevant 
correlation functions of the contribution of a Gaussian stochastic helical magnetic field. As it will be 
shown there are contributions due to the helical part of the spectrum to the correlation functions for 
all three modes, namely, scalar, vector and tensor modes. In section 4 the perturbation equations for 
the vector and tensor modes in the presence of a magnetic field are presented in the gauge-invariant 
formalism [21]. The corresponding equations for the scalar mode can be found in [19]. Moreover, 
the initial conditions are presented. In section 5 results of the numerical calculation of the angular 
power spectra determining the autocorrelation and cross correlation functions of the temperature 
(T) and polarization, that is the E-mode (E) and B-mode (B), of the CMB are presented. Due to 
the helical nature of the magnetic field apart from the autocorrelation functions of T, E and B and 
the temperature polarization E-mode cross correlation TE there are also non vanishing correlations 
between E and B as well as T and B. The latter one, in particular, is compared for a choice of 
parameters with the WMAP7 data [22]. Section 6 contains the conclusions. 
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2 Decomposition of the magnetic field contribution 



In a flat Friedmann-Robertson-Walker the lab frame is defined locally by choosing lab coordinates 
such that dt = adr and df = adx [23] . The magnetic field is treated in the lab frame in which it is 
related to the Maxwell tensor F^j, by 

Bjjx, t) = (-ijmFjmi (2-1) 

where eijm is the totally antisymmetric symbol, with €123 = 1 

The energy-momentum tensor of the electromagnetic field measured by the fundamental ob- 
server has the form of an imperfect fluid [24] 

Tap = {P + P) UaUp + pgal3 + "^Uf^aQfi) + '^a0 (2-2) 

where u" = u^^Sq is the four- velocity of the fluid and UaU°' = —1. The heat flux vanishes in 
a pure magnetic case since it is determined by the Poynting vector. The magnetic energy density 
Pb, pressure pB and anisotropic stress 7r(B) in the lab frame are given by [24], 

PB = ^ — , Pb = gPB, 7r(B) ij = -Bi{x,T)Bj{x,T) + -B^{x,T)6ij, (2.3) 

where the anisotropic stress has only non-vanishing spatial components and the vector notation 
denotes a spatial 3-vcctor. Furthermore, the term due to the Lorentz force entering the equation of 
the baryon velocity evolution and in the tight-coupling limit the photon velocity evolution which 
is derived from V^T^,^™-* = — J" and expressed in terms of quantities in the lab frame yields to 

L{x, T) = a(^JxB^ {x, t) (2.4) 

For vanishing electric field the current J is given by a J = V x B. Therefore the components of the 
Lorentz term takes the form 

i 

In the following the magnetic field contributions, namely, the energy density, the anisotropic 
stress tensor and the Lorentz term are expanded into scalar, vector and tensor harmonic functions. 

Moreover, it is used that Bi{x,T) = i?,(.x, tq) (^f^y) and = p^o (^)^ where the index refers 
to the present epoch. The energy density of the magnetic field is written in terms of the gauge 
invariant magnetic energy contrast such that 

PB = p^^ABik)Q(-''\k,x), (2.6) 

k 

where Q{k,x) denote a set of scalar harmonic functions satisfying (A + k'^)Q^^'> = (cf. e.g. 
[21]). Moreover Ab(A;) = Ab(A;,to) and similarly the dependence on tq is omitted in the following 
expressions. The magnetic anisotropic stress is determined by 

7r(,,)(f,r)=p, (2-7) 



m=0,±l,±2 



k 
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where m = denotes the scalar part and = k '^Q\ij + the vector part is determined by 

m = ±1 and Q^^"^^ ~ i^^^f^ ~^ ^iN^^ ^'^^ tensor modes are given by m = ±2 
Expanding the Lorentz term as 



X) 



(2.8) 



m=0,±l,±2 



and using that (cf., e.g., [21]) Qf'^ = —kQf\ Q-r = \kQf^ for the scalar harmonics, Q^^}^ = 

for the vector harmonics and Q^^f^ = for the tensor harmonic functions. Then the only 
non- vanishing contributions are given by the scalar and vector part 



.(0) ^(0) _ 2 1.^(0) 



L^'\k) 



The magnetic field is written as 

Bi{x,To) = J2Bi{k)Q^''\k,x) 



(2.9) 
(2.10) 

(2.11) 



which implies 



(2.12) 



A convenient representation of the scalar, vector and tensor harmonic functions in fiat space is 
given by [25, 26, 20] 



Q^^'\k,x)i = ±^(ei±ze2),e*-" 
v2 



(2.13) 
(2.14) 



(ei ± ie2)i ® (ei ± ie2)j e^^'^ . (2.15) 



The (spatial) coordinate system defined by the unit vectors ei , e.2 and es is chosen such that 63 lies 
in the direction of fc, thus 63 = k. Moreover, in the helicity basis [8] 



(2.16) 



so that • = — 1 and • = and • A; = 0. With this choice the scalar, vector and tensor 
parts of the anisotropic stress are found to be 



.^^^\k) 



2P70 



^ ^Bi{k - q)q'B,{q) [k^ - q^) - ^ 5„(fc - q)B^{q) 



ef Bi{k - q)B, {q)^ + ( e| B, {q)BS - q)k 



2 3 
3 



^(e|)^B.(A;-,-)(e|)^5,(,-). 



(2.17) 

(2.18) 
(2.19) 
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3 Helical magnetic fields 



Magnetic helicity plays an important role in the efficiency of magnetic dynamos. It provides a 
measure of the topological structure of the magnetic field, in terms of linkage and twists of its field 
lines. It is defined as an integral over the volume V by (cf., e.g. [27], [28]) 

Hm = ^ [ A-Mx (3.1) 

where A is the gauge potential and B = V x A. The expression for the magnetic helicity (3.1) 
is gauge dependent if the normal component of the magnetic field B„ does not vanish on the 
boundary of the volume. In this case magnetic helicity is a conserved quantity in the limit of 
large conductivity. Strictly speaking magnetic helicity is not defined if Bn ^ 0. However, more 
general definitions of helicity have been formulated such as a relative helicity which is manifestly 
gauge-invariant [29]. Moreover, there are many physical situations where it is not natural to assume 
that the magnetic field vanishes on the boundary such as in the case of a stellar magnetic field or 
the magnetic field inside the horizon. In [30] a different gauge-invariant definition of the magnetic 
helicity is proposed. Similarly to magnetic helicity which describes the complexity of the magnetic 
field structure there exists the concept of kinetic helicity which determine the structure of the 
velocity field v which is important in turbulence. Kinetic helicity is defined by 

Hk = J d^xv- (V (3.2) 

Therefore, sometimes in analogy to the expression of the kinetic helicity, the quantity 

Hc = ^ J (fxB ■ (v xB^ (3.3) 

is considered as a measure of the magnetic helicity which is gauge invariant, but is not an ideal 
invariant as pointed out in [27]. This form was used, e.g., in [7, 8, 9]. It describes the (electric) 
current helicity [31]. 

Assuming the magnetic field to be a homogeneous Gaussian random field with zero mean, the 
most general form of the two-point correlation function Cij in real space which is invariant under 
rotations, but not reflections is given by [32] (cf. [33, 34]) 

Cij{xi,X2) = Cij{r) = [CL{r) - Civ(r)] ^ + CN{r)dij + CAir)eijm^, (3.4) 

where 'r = xi — X2 and r = \r\. Moreover, Cl and Cn are the longitudinal and lateral correlation 
functions. The function Ca describes the asymmetric part which vanishes in the case of homoge- 
neous, isotropic random fields, which by definition arc also invariant under reflections. In the case 
of a magnetic field the asymmetric part is related to its helicity. The functions Cl and Cn are not 
independent in the case of the two-point function of a stochastic magnetic field, since V • -B = 
[32, 33, 34]. This reduces by one the number of free functions determining the two-point correlation 
function in Fourier space. Therefore, in Fourier space the correlation function for a homogeneous, 
Gaussian magnetic field reads [32] (see also, [7, 8, 9]) 

{B*{k)Bj{k')) = 5^^,Ps{k) (5ij - kk) + S^^,PAik)ieijmkm, (3.5) 
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where Psik) is the power spectrum of the symmetric part which is related to the energy density of 
the magnetic field and PAik) is the power spectrum of the asymmetric part related to the hclicity 
of the magnetic field. Moreover a hat indicates the unit vector, so that ki = ^. Following [19] the 
powers spectra are chosen to be of the form, 

/ \ 

Ps{k,kmM) = ^sf^j Wik,km) (3.6) 



k 



PA{k,km,kL) = "^Hl^-j^j W{k,km), (3.7) 

where Ab and Ah are the amplitudes and ns and ua are the spectral indices of the symmetric and 

antisymmetric parts, respectively. The spectra have to satisfy the so-called realizability condition 
|-P^(fc)| < Ps{k) which is basically a consequence of the Schwartz inequality when applied to the 
average helicity (cf. e.g. [35]). Moreover, kL is a pivot wave number and km is the upper cut-off in 
the magnetic field spectrum due to diffusion of the magnetic field energy density on small scales. 
It is assumed that this cut-off is the same for the symmetric and asymmetric parts. Furthermore, 
W{k, kfn) is a window function. The damping of the magnetic field is determined by the Alfven 
velocity and the Silk damping scale which leads to an estimate of the maximal wave number [23] 
given by (see also [36]) 

km ^ 200.694 (^-^^ Mpc"^ (3.8) 

using the values of the best fit ACDM model of WMAP7 Qb = 0.0227/1-^ and h = 0.714 [22]. The 
window function is assumed to be Gaussian of the form [19] 

Wik,km) = TT-'^ k-^e' (^) (3.9) 

where the normalization is chosen such that J d^kW{k, km) = 1. A different choice of window 
function, namely a step function, was used in [16]-[18]. 

In the continuum limit — )■ / the magnetic energy density today smoothed over the 
magnetic diffusion scale is given by 

which is valid for ns > —3, where we have used that pB = {B{x,t)^)/2. The average helicity 
measures Hm and He result in 

^« = j^r^T^("-^) (3.11) 

which is valid for > — 2 

2-K2 \kL, 

which requires ua > —4. Therefore the amplitude of the spectral function of the asymmetric part 
of the two-point function of the magnetic field can be written as 

AH = 27rlnB(^] (3.13) 



2TrV^km V kL 



Ha = ^fM"\f!^, (3.12) 
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where 



^ 1 HmKJT (^^^) magnetic helicity 
^ \ Hck^/V^"^^) current helicity 

In the numerical solutions we consider the maximal allowed contribution of the asymmetric part 
of the spectrum, that is for ua — ^5 > 0, the condition PAi^max) = Ps{kmax) is imposed, where 
kmax is the maximal wave number considered. In the opposite case, for ua — < 0, the condi- 
tion PAikmin) = Ps{kmin) is imposed, where kmin is the minimal wave number considered in the 
numerical solution. This leads to 



P70/ \Pio) r2(n^) \k 



(3.15) 



where q = kmax (kmin) for — > (< 0). Therefore, the larger the absolute value of the 
difference between the spectral indices of the asymmetric and the symmetric part of the magnetic 
field spectrum the stronger the helical contribution is suppressed. 

The two-point correlation functions of two random fields F and G in fc-space can be written in 
terms of the dimensionless spectrum Vfg as 

(PlGp) = ^VMk)S^krk'- (3-16) 

All correlation functions will be calculated in the continuum limit. The autocorrelation function 
of the magnetic energy density contrast is found to be 

[r[liS±3)Y \P^,oJ \kmj Jo 



dxe yT^) (1 _ 2zx + z^) 2 (1 + - Azx + 2z'^) 

1 

/2 / i. \ 2(nA+3) / \2 /.oo „/ ^2^2 



2/^7,0 \kmj Jo 

dxe (l-2zx + z^) 2 (x-z), (3.17) 



which reduces to the known correlation function for non helical magnetic fields (cf., e.g., [19]). 

Moreover, x = ^ and = f • Therefore, it is found that the contribution due to the asymmetric 
part of the correlation function of the magnetic field does not vanish. Similar to the case of the 
tensor modes, it is the product of two factor involving the helical part which contributes. Behind 

this is the observation that the product of two odd parity quantities results in one with even parity 
[8] . The cross correlation function between the magnetic energy density contrast and the anisotropic 
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stress in the scalar sector is given by 

1 k \^ ".9-2 

dare (l _ 2zx + z^) ^ (-1 + z'^ + zx-{l + 3z^)x^ + 3zx^) 



Jo 



c?a;e V*^^ (l - 2za; + z^) (z + 2x - Szx^) . (3. 



18) 



The autocorrelation function of TTg ^ is determined by 



'B "B 



^2 



dxe (1 _ 2zx + z^) (1 + + 2zx + (1 - 12z2)x^ - Gzx^ + 'dz^x*) 

dxe (1 _ 2za: + z^) ^ (4z + 2x - 6za;2) . (3.19) 

Since the magnetic field is helical and thus its spectral function has an asymmetric part there 
are two different, relevant two-point correlation functions for the anisotropic stress of vector and 
tensor modes. The symmetric part which determines the angular power spectra of the CMB due 
to all parity even or all parity odd modes is given for the vector modes by 



72 



^y/^y<--) (^) 



[r(ii^)]^ \pjoJ \h 

J dzz"=+2e"K^) J dxe^^) (1 - 2zx + z^) (1 - x^)(\ + - ?,zx + Iz^x^) 

P-1O \kmj Jo 



j dxe'^^-^^ ^'^{l — 2zx 



+ Z^)2{l-X^) 



(3.20) 



This correlation function is used to calculate the angular power spectra of the temperature and 
polarization autocorrelations, that is Cf^, Cf^ and Cf^, and the temperature polarization E- 

mode cross correlation Cj^ . The two-point correlation function for the anisotropic stress of vector 
modes determining the angular power spectra of the CMB due to a mixture of parity even and 
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parity odd modes is found to be 



27r2 



36 



6+ns+nA 



xe 



-(^y 

\km J 



+2 



dxe 



'''{\-2zx + z^)^{l-2zx){l-x'^) 



-f 

Jo 



dzz'^^+^e 



dxe 



\l - 2zx + z^) 2 (1 - 2zx){l - x^) 



(3.21) 



which is used to calculate the parity-odd correlations of the CMB anisotropics, that is Cf^ and 
Cj^ . It is interesting to note that in the case = n^ + l the correlation function (3.21) identically 
vanishes. 

For the tensor modes the two-point correlation functions of two even-parity modes is given by 



(7r(+2)*(^)4+^)(fcO + vri-^)*(^)4"^^(fc')) = ~—k^, 



24 



( ( A 



r2(n^) \p^^) \k, 



2(3+n5) 



\km) 



X / dzz'^s+^e 

10 

\PjoJ \k, 



f 

Jo 

fnB\ 

P-io) 

/"°° dzz"^+2e"K^) 
Jo 



dxe 



/: 



\l-2zx + z'^)2{l + x^) 



l-2zx + —{1 + x^) 



k 



2(3+nA) 



dxe 



2(^) ^^^1 o.^ , .2^:iA.-l 



'{l — 2zx + z^) 2 {l — zx)x 



(3.22) 



which determines Cj"^, Cf^, Cf^ and C|"^ for tensor modes. The parity-odd correlations of the 
CMB anisotropies due to tensor perturbations arc determined by 



2tt^ 

fe3 



xe 



12 



r(^) \PioJ Kp^oJ \km 



k 



\ 6+ns+nA 



■J 



\kmj 



Jo 



dzz 



-(^y 

\km) 



+ 



f 

Jo 



dzz'^s+'e 



dxe 



+2 -2U!r ^= 



'^x(l - 2zx + z^) 2 



l-2zx + —{l + x^) 



dxe^^''"^^ 



(1 - 2zx + z^)~^{\ - zx){\^x^ 



(3.23) 



yielding to non-vanishing cross correlations between the E- and B- mode Cf^ and the temperature 
and B-mode Cj^ . These expressions agree with [8]. The dimensionless spectra (cf. equation 
(3.16)) determining the correlation functions for different choices of ns and ua arc shown in figure 
1 for the scalar modes and in figure 2 for the even-parity correlations of vector and tensor modes. 
Finally in figure 3 the odd-parity correlations of vector and tensor modes are reported. In the 
numerical solutions the amplitude of the helical component is taken to be the maximal allowed 
value allowed by the rcalizability condition (cf. equation (3.15)). As can be seen from figure 1, 
as in the case of nonhelical magnetic fields the magnetic energy density and anisotropic stress are 
anticorrelated in the case of scalar modes. Moreover, for unequal spectral indices for the symmetric 
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Figure 1: Left: The spectra determining the correlation functions for scalar perturbations for 
B = 10 nG, ns = —2.9 and ha = —1.9 and for comparison the case of a nonhelical magnetic field 
has been included for B = 10 nG, ns = —2.9. Right: The spectra determining the correlation 
functions for scalar perturbations for B = 10 nG, ns = —2.9 and n^ = —2.9. The amplitude of the 
helical component is taken to be the maximal allowed value allowed by the realizability condition. 
Moreover, the maximal value of the wave number is set to kmax/km = 100. 



and asymmetric parts of the magnetic field correlation function there is a strong suppression of 
the helical contribution for the scalar modes. This manifests itself in the small change between the 
solutions for n^ = —1.9 and a nonhelical magnetic field in figure 1 (left). As can be appreciated 
from figure 2 for n^ ns the helical part is strongly suppressed in the even parity correlation 
functions for vector and tensor modes which is also observed in the case of scalar modes. The non- 
vanishing odd-parity correlations for the vector and tensor modes in figure 3 constitute a distinctive 
feature of the helical nature of the magnetic field which lead to the non- vanishing cross correlations 
Cf^ and Cj^ in the CMB. 

4 Perturbation equations 

In this section the perturbation equations and the initial conditions in the presence of a magnetic 
field for the vector and tensor modes in the gauge-invariant formalism are presented. Perturbations 
are considered around a flat background ds^ = a?{T){—dT'^ + 6ijdx^dx^), where o(r) is the scale 
factor. The corresponding equations for the scalar sector can be found in [19]. However, we 
will comment briefly on the initial conditions. Current numerical codes to calculate the CMB 
anisotropies, such as [11, 12, 13], set the initial long after neutrino decoupling at which time there 
is already a non vanishing value of the neutrino anisotropic stress. However, in most models 
the magnetic field will be generated long before neutrino decoupling. This leads to the magnetic 
field anisotropic stress being the only contribution to the total anisotropic stress. As shown in 
[37, 38, 18] in the case of scalar modes the magnetic anisotropic stress causes the comoving curvature 
perturbation on superhorizon scales to evolve approaching a constant value shortly after neutrino 
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Figure 2: Left: The spectra determining the even parity correlation functions for vector and tensor 
modes for B = W nG, Ug = —2.9 for ua = —1.9 and for comparison the case of a nonhehcal 
magnetic field has been included for B = W nG, ns = —2.9. Right: The spectra determining 
the even parity correlation functions for vector and tensor modes for B = 10 nG, ns = —2.9 and 
riA = —2.9. The amplitude of the helical component is taken to be the maximal allowed value 
allowed by the realizability condition. Moreover, the maximal value of the wave number is set to 

^max/^m — 100. 
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Figure 3: The spectra determining the odd-parity correlation functions for vector and tensor modes 
for i? = 10 nG, ng = —2.9 and ha = —1.9 {left) and ua = —2.9 [right). The amplitude of the 
helical component is taken to be the maximal allowed value allowed by the realizability condition. 
Moreover, the maximal value of the wave number is set to kmax/km = 100. 
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coupling. Moreover, at this time the neutrino anisotropic stress approaches the compensating 
solution, in which it balances the magnetic anisotropic stress [39, 16]. However, the comoving 
curvature perturbation has acquired an additional contribution due to its evolution before neutrino 
decoupling. This is in addition to any primordial curvature perturbation generated, for example, 
during inflation. The additional contribution is determined by [18] 

Cc.-\lt,nPln^, (4.1) 

6 TB 

where = ^ and tb is the time of (instantaneous) generation of the magnetic field, which 
would correspond to the time of the phase transition (e.g. [5]) or to reheating if generated during 
inflation (e.g. [4]). In [18] a distinction was made between the compensated magnetic mode with 
= and the adiabatic-like passive mode with ( given by equation (4.1). In [38] this distinction 
was not made and here this approach is followed. The initial conditions used in the numerical 
solution include as a new parameter P = ln^. 

4.1 Vector perturbations 

There are two gauge-invariant equations [21]. The amplitude of the shear of the normal vector field 
to the constant time hypersurface, ag^\ is determined by 

+ 2WaW = k [n, (ttW + ttW) + , (4.2) 

where a dot indicates the derivative with respect to conformal time r and H = ^. The amplitude 
of the vorticity of the matter velocity field V^^^ is gauge-invariant and satisfies [21] 

+ (1 _ 3c2)Wy(i) = -^-^7r(i), (4.3) 

where w determines the equation of state of the fluid P = wp with P the pressure and p the energy 
density and vr^^^ denotes the anisotropic stress. Moreover, Cg is the adiabatic sound speed. For 
massless neutrinos (v) and cold dark matter (c) equation (4.3) yields to 

K« = ~4'^ (4.4) 
l/i^^+^FW = 0. (4.5) 

Deep inside the radiation dominated era photons and electrons arc tightly coupled through Thomson 
scattering and the latter are tightly coupled with the baryons via Coulomb interaction. Thus the 
baryon, electron and photon fluids are well described by a single fluid. In the tight coupling limit 
the vorticity fields of photons (7) and baryons (b) are determined by 

VJ}'^ = (4.6) 
(1) _ 



'{y^''-yi'')-^-'S. (4.7) 



where ^ is the mean free path of photons between scatterings which is determined by the number 
density of free electrons ng and the Thomson cross section ctt, = angUT- Moreover, in the 
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baryon vorticity equation (4.7) R = |^ and the magnetic field contribution is due to the vector 
component of the Lorentz force (cf. equation (2.10)). At very early times the mean free time 
between scatterings is much smaller than the Hubble time implying a comparatively large value of 
T~^. This leads to problems in the numerical integration and was first solved in the case of scalar 
perturbations by using an iterative solution at early times and the exact equations at later times 
[40, 41, 42]. For the vector perturbations we use a similar approach which results in 



K 



vr + 

R k 



R 



~ l + i^s""^ 1 + RV ^^^b J' 



(4.8) 
(4.9) 



where the shift V^^^ = V^^"^^ — Vy^^ is determined by 



V(i) = 



1 + 2.^-^ 



H -1 



l + R 

+ ?(v;'"-K,w) 



l + R 



(4.10) 



and in the tight-coupling limit the term V^^^ — V^''' is neglected. 

As pointed out for the scalar modes, initial conditions for the numerical solutions are usually set 
after neutrino decoupling when the neutrino anisotropic stress is already non vanishing. Therefore 
before neutrino decoupling the magnetic field is the only source of anisotropic stress [37, 38, 18]. 
In the case of the scalar perturbations this leads to an additional contribution to the comoving 
curvature perturbation as discussed before. Thus the only source for the amplitude of the shear of 
the normal vector field to the constant time hypersurfaces ag^^ is given by the anisotropic stress of 
the magnetic field, therefore equation (4.2) implies that at Ti, 



a, 



(1) 



(r.) 



(7, 



(1) 



irB) + 



(1) 



1 - 



(4.11) 



where tb is the time of generation of the magnetic field. In case it is generated during inflation this 
time is chosen to coincide with the beginning of the radiation dominated era (see also section 4.2). 
After neutrino decoupling, ignoring any photon anisotropic stress, and neglecting the contribution 
from the multipole £ = 4 the relevant equations from the Boltzmann hierarchy in addition to 
equation (4.4) are given by [20] 



TT, 



(1)' 



35 



(1)' 



(1) 



24 



(4.12) 
(4.13) 



where A^2^^ = ^^'^i'^^- This gives together with equation (4.2) 

^(1)"' + 3^(1)" 1 /99 8J^\ /9 _ 8a.\ 

^ ^ X " ^5 V35 ^ x-' ) 5x2 ) 



X 



(1) 



8 l^/yTT^ 



(4.14) 
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where a prime denotes the derivative with respect to x = kr. On superhorizon scales a; <C 1 
equation (4.14) is solved by 



4'\r) = -':^^'^> + {kr) 



Cie-2 



y/32n„/5-lln(kT) 



+ 0262 



§-^320^/5-1 In(fcT) 



(4.15) 



so that for fixed k and late times the compensating solution is approached, tt, 

implies the vanishing of the righthandside of equation (4.2) leading to ag^^ ^ and therefore 
suppressing any contribution due to the evolution before neutrino decoupling. Thus, finally, the 

compensating magnetic initial conditions for the neutrino anisotropic stress together with equations 
(4.2), (4.4), (4.5), (4.8), (4.9), (4.12) and (4.13) yield the following overah initial conditions in terms 
of a; = kr, 



a, 



(1) 



15 n^TT^j^^ 

14 15 + 4f], 



■X, 



V, 



(1) 



^1 



IT 



(1) 



-X, 



(1) 



7^(1) 



TT 



-1 + 



45 



14 Ib + AQy 



(1) 



T/(l) - -^J^K 

"8 0/«^' 

(1) 



TT, 



(1) 







0,^ TTj^ 



-1 + 



15 



14 15 + 40.^ 



(4.16) 



which agree with those of [18]. In deriving these initial conditions it was used that for the numerical 
solutions to calculate the CMB anisotropics these are set deep inside the radiation dominated era 
when the baryon-photon fluid is tightly coupled. 



4.2 Tensor perturbations 

(2) 

The metric tensor perturbations are determined by one gauge-invariant amplitude Hj, which 
satisfies [21] 



(2) 
7 



(4.17) 



It is useful to recall that the initial conditions for the numerical solutions of the Boltzmann equations 
are set after neutrino decoupling, that is at an initial time Ti > t^. Therefore, as in the case of the 

(T) 

scalar and the vector perturbations, the evolution of Hj, ' has to be determined before neutrino 
decoupling, when the magnetic anisotropic stress is the only source in equation (4.17), and after 
neutrino decoupling when also the anisotropic stress of the neutrinos contributes. On superhorizon 

(2) 

scales, X = kr <^ 1 using the evolution of Hr^ during the era before neutrino decoupling, at a 

(2) 

regular solution for Hj, is found to be 



h: 



(2) 



(r,)~iji.')(rB) + Q^7rg) In^, 



TB 



(4.18) 



where Tb is the time of generation of the primordial magnetic field if this takes place after the 



0. If, 



(2) 

beginning of the radiation dominated era. In this case it is natural to assume Hj. (tb) 

(2) 

however, the magnetic field is generated during inflation, then the evolution of Hj^ during inflation 
has to be matched to the evolution during the radiation dominated era before neutrino decoupling 
which is considered here. However, this is beyond the scope of this article and will be considered 
elsewhere (for work in this direction for the scalar modes see [43]). For simplicity, it is assumed 
that Tb = trh the time of reheating when the standard radiation dominated era begins and in the 
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(2) (2) 

numerical solution Hj. (tb) = 0. After neutrino decoupling the amplitude of the tensor mode Hj, ' 
and the neutrino anisotropic stress ir^^ satisfy [20], 

2 



H, 



(2)" 



+ -i?^' + H, 



r(2) 



TT 



(2)' _ 



at; 



(2)' 



X 



^ '1 jy 1 



(2) 



7^5 



(2) 



(4.19) 
(4.20) 

(4.21) 



where a prime indicates the derivative with respect to a; = A;t and the neutrino multipole for £ = 4 
has been neglected. Equations (4.19) to (4.21) can be combined to give a fourth order differential 
equation for h'^\ 



(2)"" 6 (2)'- 



+ 



8 2 15 + AQ.^ 

"11 ~i~ "3" n 



H. 



(2)' 



, 30 (2)' 



1 2 



(2) 



1 Vt^TT%^ 



(4.22) 



On superhorizon scales x <C 1 a rcgTilar solution at some time > r^, which is going to be taken 
the initial time to start the numerical evolution of the Boltzmann code, is given by 



<)(r,) 



1 



hx 



2(15 + 4fJ^)_ 
5x^ 



+ Q-yVr^^ In — 



TB 



1 



5.T 



+ 0{x'), 



(4.23) 



28(15 + 4^]^) 

where the solution has been determined by matching the solution before and after neutrino decou- 
pling at r = Tn- Similarly, using equations (4.19) to (4.21) the evolution of the anisotropic stress 
of the neutrinos is determined by. 



1-4 ri=> 



^11 



+ 



+ 



X 

2 1 
35 X 
1 

7 ~ 



4 ) 

x^ 



+ 



1 



x'- 



44 m, 



7 



16 ill, 

y x4 



TT. 



(2)" 



4(5 + 280^) 1120 



6 48fi^ 
7x2^ 5x4 



+ 



x^ 



TT, 



(2)' 



32fi,. 



5x6 



^(2) = 

5 



x^ 



(4.24) 



On large scales x <C 1 a regular solution is found to be 



4^^(r.) = - 



__^7^(2) 



+ 



(2) 



15 + 40^ ^ 



.(2) 



In^ + 

15 + TB 14 VLy 15 + 40^ 



x2 + 0(x^)(4.25) 



where equations (4.20) and (4.21) have been used to determine the free constant. Therefore, the 
solution approaches the one corresponding to the so-called compensating mode. However, since 
there are no independent constants multiplying the passive and the compensating mode the two 
parts of the initial conditions are not treated separately as in [18]. This follows [38] where in the case 
of scalar perturbations no separation into a passive and a compensating mode was made. We close 
this section by noting that in the numerical solution the standard tight-coupling approximation is 



used [18], in the notation of [20], 6 



(2) 



TT- 



(2) 



photons, and = —^Q^\ b'^^ = for the polarization. 



^tcH!^^ for the anisotropic stress of the 
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5 Results 



In [19] the CMB temperature anisotropies and polarization due to scalar perturbations in the 
gauge-invariant formalism in the presence of a stochastic magnetic field have been calculated using 
a modified version of CMBEASY [42]. As opposed to the case considered here, it was assumed 
that the magnetic field is non helical. The calculation of the angular power spectra of the CMB 
temperature anisotropies and polarization in the presence of a helical stochastic magnetic field as 
described by the two-point correlation function (3.5) has been done by expanding the numerical 
code of [19]. In the case of the scalar perturbations the initial conditions have been changed 
in order to include the contribution of the magnetic field to the curvature perturbation due to its 
evolution before neutrino decoupling. Moreover, a new part has been added to the modified version 
of CMBEASY [19] to include the numerical solution of the Boltzmann hierarchy and the calculation 
of the CMB anisotropies and polarization for the vector and tensor modes using the total angular 
momentum approach of Hu and White [20] . In the numerical solution here it is not assumed that 
there is an explicit separation into a magnetic mode and a passive mode as in [18]. The initial 
conditions for the numerical solution are set long after neutrino decoupling. For scalar and tensor 
modes a new parameter, /3 = In is included encoding the evolution of the comoving curvature 
perturbation in the case of the scalar mode and the amplitude of the tensor mode, respectively, 
after the creation of the magnetic field at tb and before neutrino decoupling at t,^. The vector 
mode is not affected significantly by the presence of a magnetic field before neutrino decoupling. 
The spectrum of the stochastic magnetic field is effectively cut off at the magnetic diffusion scale 
using a Gaussian window function [19]. The contribution of the magnetic field is determined by the 
two-point correlation functions involving the magnetic energy density and anisotropic stress leading 
to convolution integrals in Fourier space. Rather than using an approximation, as in previous work, 
which however, also differs in the sharp cut-off of the magnetic field spectrum as opposed to the 
Gaussian window function used here, e.g. [16, 17, 18], these integrals are calculated numerically in 
the code as was done in [19] for the non helical case for the scalar mode. The angular power spectra 
describing the CMB temperature anisotropies and polarization are obtained following a treatment 
similar to the one of correlated isocurvature modes (e.g. [44]) [18, 19] with the relevant correlation 
functions being the ones of the magnetic field contribution. 

In figures 4-6 the angular power spectra determining the autocorrelation and cross correlation 
functions of the temperature anisotropies and polarization of the CMB due to the scalar, vector 
and tensor modes are shown for a choice of the magnetic field parameters and the best-fit values 
of the 6-parameter ACDM model of WMAP7 [22], in particular, = 0.0445, = 0.738 and the 
reionization optical depth r = 0.086. In all numerical solutions the amplitude of the helical part 
of the magnetic field is taken to be the maximally allowed value by the realizability condition (cf. 
equation (3.15)). Moreover, for ha = —1.9 the maximal wave number is set to kmax = lO^^m- 
The spectral indices of the symmetric and asymmetric parts are chosen to be negative. Whereas 
magnetic fields generated during infiation with cosmologically relevant field strengths result to 
have negative spectral indices (e.g. [4]) those generated during a phase transition, such as e.g. the 
electroweak one, have to satisfy causality constraints which require positive spectral indices [45, 8]. 
The asymmetric part of magnetic field correlation function causes non vanishing cross correlations 
between the E-mode and B-mode and the temperature and the B-mode, respectively. These are 
shown in figure 5. Current observations of the B-mode of polarization are consistent with zero 
[22, 46]. Thus comparing in particular the cross correlation between temperature and the B-mode 
seems a promising possibility to constrain the helical contribution. In figure 6 the TB- and the BB- 
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B=5 nG, ns=-2.9, |3=0 B=5 nG, ns=-2.9, |3=0 




B=5 nG, n=.=-2.9, p=0 

B=5 nG, ns=-2.9, p=0 




Figure 4: The TT, EE, TE and BB angular power spectra for the scalar, vector and tensor modes 
for the magnetic field strength B = 5 nG, spectral index of the symmetric part of the magnetic 
field correlation function ns = —2.9. The spectral index of the asymmetric part is assumed to be 
UA = —2.9 and ua = —1.9, respectively. The pure magnetic mode is shown, /3 = 0. 

angular spectra have been plotted for different choices of parameters together with observational 
data from WMAP7 [47]. As can be appreciated from figure 6 the constraint on the contribution to 
Cf^ due to the vector mode is at high values of the multipoles i and for the tensor modes at low 
values of £. Whereas the vector modes do not depend on the parameter /3 encoding the evolution 
of relevant quantities upto the time of neutrino decoupling, a nonzero value of /3 constrains the 
magnetic field strength to lower values due to the tensor modes. Assuming that the magnetic 
field is created during inflation and thus setting tb to the beginning of the standard radiation 
dominated era, so that /3 = In then assuming reheating at 10^*^ GeV and together with Ty = 1 
MeV resulting in (3 = 30, puts an upper limit on the magnetic field strength oi B = 2 nG for 
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Figure 5: The angular power spectra Cf" and due to vector and tensor modes for the 

magnetic field strength B = 5 nG, spectral index of the symmetric part of the magnetic field 
correlation function ns = —2.9. The spectral index of the asymmetric part is assumed to be 
UA = —2.9 and = —1.9, respectively. The pure magnetic mode is shown, /3 = 0. 
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Figure 6: The TB and BB angular power spectra for different parameters comparing with WMAP7 
data. Note that for TB the plotted angular power spectrum is {i + l)Cj^/(27r) in accordance with 
the WMAP7 data and not £{£ + 1)CJ^/{2t:). 



IT'S = ''T'A = —2.9. For lower reheat temperatures larger values of the magnetic field strength are 
allowed. The WMAP7 data for the BB spectrum are less constraining as can be seen in figure 6 
{right). In the numerical solutions it was assumed that there is no correlation between the magnetic 
field contributions and any primordial curvature perturbation or tensor modes from, e.g., inflation. 



18 



However, if the magnetic field is generated during inflation then one might expect a correlation 
which deserves further study. This has recently been considered in [48]. 

6 Conclusions 

The CMB anisotropies and polarization in the presence of a stochastic helical magnetic field have 
been calculated for scalar, vector and tensor modes. For this purpose the modified version of CM- 
BEASY [19] calculating the CMB anisotropies due to the scalar perturbations in the presence of a 
non helical magnetic field has been expanded. Firstly the numerical solution for the corresponding 
correlation functions has been included for scalar, vector and tensor modes. Secondly a new part 
has been added to include the calculation of the CMB anisotropies and polarization due to vector 
and tensor modes using the total angular momentum approach of Hu and White [20] . A Gaussian 
window function is used to effectively cut-off the magnetic field spectrum at a wave number cor- 
responding to the magnetic damping scale. In the case of the scalar and tensor perturbations the 
initial conditions for the numerical solution which are set long after neutrino decoupling include 
a contribution encoding the evolution of relevant quantities due to the presence of the magnetic 
field before neutrino decoupling. In the case of the scalar perturbations the comoving curvature 
perturbation grows upto the time of neutrino decoupling due the magnetic anisotropic stress. Af- 
ter neutrino decoupling the neutrino anisotropic stress compensates the magnetic anisotropic stress 
and the comoving curvature perturbation becomes a constant on superhorizon scales [37, 18, 38]. 
In the case of the tensor modes a similar behaviour is found for the amplitude of the tensor mode 
which has been shown here explicitly. The presence of a magnetic field prior to neutrino decou- 
pling does not affect significantly the evolution of the vector modes. Using a standard Boltzmann 
solver the contribution due to presence of the magnetic field before neutrino decoupling has been 
included in the initial conditions for the numerical solution which are set long after neutrino de- 
coupling. However, since the presence of a magnetic field affects the evolution of the scalar and 
tensor perturbations before neutrino decoupling in order to be more precise one would have to start 
the numerical evolution and thus set the initial conditions for the Boltzmann code before neutrino 
decoupling. A problem we hope to address in the future. 

In the case of a helical magnetic field in addition to the temperature (T) and polarization E- 
and B- mode autocorrelation spectra and cross correlation TE angular power spectrum there are 
also the cross correlation EB and TB angular power spectra. The latter one has been used to 
illustrate the use of the current WMAP7 data to constrain the magnetic field parameters. 

In comparison to earlier work on the effects of a helical magnetic field [7, 8, 9] on the CMB here 
a full numerical treatment has been provided including the correct initial conditions and evolution 
equations in the presence of a magnetic field as well as including numerical solutions for the different 
correlation functions of the magnetic field contributions. Moreover, the magnetic field spectrum is 
effectively cut-off using a Gaussian window function. 
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